Two-photon Indirect Optical Injection and Two-color Coherent Control in Bulk 
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Using an empirical pseudopotential description of electron states and an adiabatic bond charge 
model for phonon states in bulk silicon, we theoretically investigate two-photon indirect optical 
injection of carriers and spins and two-color coherent control of the motion of the injected carriers 
and spins. For two-photon indirect carrier and spin injection, we identify the selection rules of band 
edge transitions, the injection in each conduction band valley, and the injection from each phonon 
branch at 4 K and 300 K. At 4 K, the TA phonon-assisted transitions dominate the injection at 
low photon energies, and the TO phonon-assisted at high photon energies. At 300 K, the former 
dominates at all photon energies of interest. The carrier injection shows anisotropy and linear- 
circular dichroism with respect to light propagation direction. For light propagating along the (001) 
direction, the carrier injection exhibits valley anisotropy, and the injection into the Z conduction 
band valley is larger than that into the X/Y valleys. For a~ light propagating along the (001) ((111)) 
direction, the degree of spin polarization gives a maximum value about 20% (6%) at 4 K and —10% 
(20%) at 300 K, and at both temperature shows abundant structure near the injection edges due to 
contributions from different phonon branches. For the two-color coherent current injection with an 
incident optical field composed of a fundamental frequency and its second harmonic, the response 
tensors of the electron (hole) charge and spin currents are calculated at 4 K and 300 K. We show 
the current control for three different polarization scenarios: For co-circularly polarized beams, the 
direction of the charge current and the polarization direction of the spin current can be controlled 
by a relative-phase parameter; for the co-linearly and cross- linearly polarized beams, the current 
amplitude can be controlled by that parameter. The spectral dependence of the maximum swarm 
velocity shows that the direction of charge current reverses under increase in photon energy. 

PACS numbers: 42.65.-k,72.25.Fe,72.20. Jv,78.20.-e 



I. INTRODUCTION 

Silicon is a dominant material in the microelec- 
tronics industry. It has also attracted much atten- 
tion in optoelectronics^— due to its low absorption at 
telecommunication wavelengths near 1.55 /xm, and in 
spintronics^— due to its long spin transport length and 
spin relaxation time£~— In both fields, a full under- 
standing of the optical properties in bulk silicon is very 
important for further applications. Optical methods 
can provide an effective way to generate carriers and 
spins in semiconductors, to control^— their motions 
by the phase coherence of different components of in- 
cident laser beams, and to detect the properties of car- 
riers and spins H Recently, the direct detection of spin 
currents using second-order nonlinear optical effects has 
been proposed^ and realized experimentally^ 

Because silicon is an indirect gap semiconductor, with 
an indirect gap Ei g = 1.17 eV and a direct gap E g — 
3.4 eV^Si there is a degenerate indirect "£" -photon opti- 
cal transition assisted by phonon emission or absorption 
at £hu! < Eg. This optical response is about two orders 
of magnitude weaker than that in direct gap semiconduc- 
tors. While the weak response results in low loss, which 
is important in realizing optoelectronics devices, it can 
make optical coherent control less effective. 

By using circularly polarized light, spin polarized car- 



riers can be injected^ Generally, one- and two-photon 
injection are the most widely used schemes. For coherent 
current control, the minimum requirements depend on 
the semiconductor crystal structures: For low symmetry 
semiconductor structures with nonvanishing second or- 
der nonlinearity, such as the wurtzite structure^ current 
can be injected by even a single frequency laser beam; for 
high symmetry semiconductor structures with vanishing 
second order nonlinearity but non-vanishing third-order 
nonlinearity, such as the diamond structures, current in- 
jection requires at least a two-color laser pulse with one 
fundamental frequency and its £ th harmonic ("1 + £" ef- 
fects). The control parameters are taken as a relative- 
phase parameter between Cartesian components or be- 
tween the frequency components of the two-color laser 
beams. However, most coherent control studies to date 
have focused on absorption across the direct gapf 2 ^— 
even when considering the indirect gap semiconductors 25 ; 
seldom has coherent control by absorption across an in- 
direct gap been considered* 2 ^— due to the weak optical 
response. For silicon, which has diamond structure and 
vanishing second-order nonlinearity, it is only the second 
of the coherent control schemes mentioned above that is 
applicable. 

For two-photon indirect optical carrier injection in 
bulk silicon, most experimental studies have focused 
on the two-photon absorption coefficient 2 ^— and its 
anisotropy, 33 which is important in optoelectronics de- 
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vices; theoretical studies 3 ^— are mostly based on the 
parabolic band approximation and on a phenomenologi- 
cal electron-phonon interaction. For the current injection 
by coherent control, Costa et al£l and Spasenovic et al2^ 
used THz radiation to detect "1+2" injected current in 
bulk silicon, and confirmed that the current can be con- 
trolled by the phase parameter of the laser beams. Zhao 
and Smir]™. measured the time- and space-evolution of 
the indirect optical injected electrons and holes by phase- 
dependent differential transmission techniques. Yet a full 
band structure investigations of the two-photon indirect 
optical injection of spins and spin current are still lacking. 

Previously we studied the one-photon indirect optical 
injection of carriers and spins ^ and the spectral de- 
pendence of the two-photon indirect absorption coeffi- 
cients and their phonon-resolved injection rates at 4 K 
and 300 In this paper, we continue the study of 

the two-photon indirect optical injection of carriers and 
spins, and consider as well the coherent control of the 
injected charge and spin currents by "1+2" effects. We 
present the detailed results of two-photon indirect carrier 
and spin injection under o~ light propagating along (001) 
and (111) directions; due to the symmetries of bulk sili- 
con, the injection with cr + light has the same carrier and 
spin injection as with a~ light, but with the opposite spin 
polarization. The injection in each conduction band val- 
ley, the anisotropy and the linear-circular dichroism with 
respect to light propagation direction, the corresponding 
phonon-resolved spectra, and the degree of spin polar- 
ization (DSP) are discussed. We also consider the coher- 
ent control of the motion of optically injected electrons 
and holes under particular two-color optical fields: co- 
circularly polarized beams, co-linearly polarized beams, 
and cross-linearly polarized beams. 

In optical absorption, the electron-hole interaction 
plays an important role especially in determining the cor- 
rect absorption edges. First principle studies^ of the di- 
rect gap optical absorption shows that the excitonic effect 
can strongly change the lineshape even for high photon 
energies in silicon. For indirect one— and two-photon 
injection ) 37 i 38 investigations within the parabolic band 
approximation show that this neglect does not change 
the absorption lineshapes at energies more than a few 
binding energies above the band gap; however, a full 
band structure investigation is still lacking due to dif- 
ficulty in numerical calculation of the wave functions of 
the electron- hole pair. In this paper, as a preliminary 
calculation, we neglect the excitonic effect. 

We organize the paper as follows. Two-photon indi- 
rect carrier and spin injection are presented in Sec. [TXJ 
In this section, we first describe a perturbation model 
for two-photon indirect optical injection, and then give 
the numerical results under an empirical pseudopotential 
model for electronic states and an adiabatic bond charge 
model for phonon states. In Sec. IIII1 we study the in- 
terference current injection under a two-color laser beam 
and the coherent control. We conclude in Sec. IIVI 



II. TWO-PHOTON INDIRECT CARRIER AND 
SPIN INJECTION 

A. Model for two-photon indirect injection 

For an incident laser beam with electric field E(t) = 
E^e~ lult + c.c, the two-photon optical injection rates of 
electrons and their spins are generally written as 
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From these rates, the actual injected carrier density and 
spin density can be calculated once the pulse duration 
is specified. In this paper, superscripts indicate Carte- 
sian coordinates, and repeated superscripts are to be 
summed over. For bulk silicon, the lowest conduction 
band has six equivalent valleys, which are usually de- 
noted as X, X, Y, Y, Z, Z. The two-photon indirect tran- 
sitions have the same initial and final states as that of 
one-photon indirect transitions^ The injection coeffi- 
cients can be written as the form j{ abcd = J^j Aj bcd with 
j\abcd identifying the injection into the I th valley. Fermi's 
golden rule gives A abcd — J2 CV \± ^Vcv\± w ^ n 
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The coefficient Aj b ^x± gives the contribution to the in- 
jection by indirect optical transition between conduction 
band c and valence band v, with the assistance of an emit- 
ted (+) or absorbed (— ) phonon in the X th -mode; there 
are two modes each for the transverse acoustic (TA) and 
optical (TO) branches, and one mode each for the lon- 
gitudinal acoustic (LA) and optical (LO) branches. The 
operator A in Eq. stands for the identity operator in 
carrier injection, and the f th component of spin operator 
in spin injection. The optical transition matrix elements 
are given as 
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where e 



|. Here k c and k v are the electron and hole 
wave vectors respectively; c = {c, <r c }, c' = {c, <j' c }, and 
v = {v, a v } are full band indexes with a c , a' c and a v being 
the spin indexes; n and rh are band indices for interme- 
diate states; \ck c ) and e c k are the electron states and its 
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energy respectively; and uj nm (k) is denned as huj nm (k) = 
£jik — £jnk- The phonon energy is given by hfl q \ for 
wavevector q and mode A, the equilibrium phonon num- 
ber is N q \, and N qX ± = N qX + ^±^. The velocity matrix 
elements are %„-,(fc) = (nk\v\fhk) with the velocity op- 
erator v = dH e /dp, and H e is the unperturbed electron 
Hamiltonian. The electron-phonon interaction is written 
as H ep = J2 q x H e x p {q){a qX + a t _ qX ) with a q \ standing for 
the phonon annihilation operator. Its matrix elements 
are M n k c mk v X = {nk c \H e x p (k c - k v )\fhk v ). 

The injection coefficient £f bcd is a fourth-order tensor 
and ^j abcd i s a fifth-order pseudotensor. Both of them 
are symmetric on exchange of indices a and b, and on 
exchange of indices c and d. They have the proper- 

ties (£f cd )* = tf dab and ^/ afccrf )* = ({ cdab . Further- 
more, time-reversal symmetry gives ^ bcd = (£j fccd )* and 
^fabed _ _ ^fabcd\ j n faj^ silicon, each conduction 

band valley has C± v symmetry. Under this symmetry, 
^abed jjg^g g j x nonzero independent components, 
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^fabed a j go j^ ag g j x nonzero independent components, 
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The injection coefficients A abcd can be obtained by a 
proper rotation operation that transforms the Z valley 
to the I th valley. Using inversion and time-reversal sym- 
metries, all £ abcd are identified as real numbers, and all 
^fabed are p Ure imaginary numbers; A] b ^ d x± shares the 

same symmetry properties as A] bcd , while A abcd belongs 
to the higher symmetry group Oh , and has fewer nonzero 
independent components 
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All these components are related to the nonzero injection 
coefficients in the Z valley by 
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With all these coefficients, the injection rates for laser 
pulse with any polarization and propagating directions 
can be evaluated. In Appendix ^ we give in detail the 
carrier and spin injection rates for circularly-polarized 
light with any propagating direction, and the carrier in- 
jection rates for linearly-polarized light with any polar- 
ization and propagating directions. In the following, we 
focus on light propagating along (001) and (111) direc- 
tions. 



B. Results 

For quantitative calculations of the two-photon indi- 
rect injection rates, a full band-structure description of 
the electron and phonon states is necessary. Here we 
use an empirical pseudopotential modeler— for electron 
states and an adiabatic bond charge model^ for phonon 
states. All the parameters used in the empirical pseu- 
dopotential model and the adiabatic bond charge model 
are the same as those in the calculation of one-photon op- 
tical spin injection^ From the empirical pseudopotential 
model, the calculated direct band gap is E g = 3.43 eV, 
the indirect band gap is Ei g ~ 1.17 eV; the band edge 

for the conduction band is located at fc° = 0.85 and 
for the valence bands at the T point, fe° = 0. From the 
adiabatic bond charge model, the energies for phonons 
with wavevector fc° are 19 (TA), 43 (LA), 53 (LO), and 
57 (TO) meV. Within the pseudopotential scheme we de- 
termine the electron-phonon interaction, and then eval- 
uate the matrix elements H x v {q) using the calculated 
electron and phonon wavefunctions. With all these quan- 
tities calculated, the two-photon indirect gap transition 
matrix elements in Eq. ^ are calculated using the lowest 
30 electron bands as intermediate states to ensure con- 
vergence. The injection coefficients given in Eq. are 
evaluated using an improved linear analytic tetrahedral 
method^ 

In our calculation, the valence bands include heavy 
hole (HH), light hole (LH), and spin split-off (SO) bands; 
the conduction bands include the lowest two conduction 
bands. Our results are shown in Fig. [1] for the spectra of 
nonzero components of ^ bcd and in Fig.[5]for the spectra 
of nonzero components ofC^ abcd at 4 K and 300 K, respec- 
tively. The full two-photon indirect gap injection rates in 
Eq. (TTJ) can be identified for any polarization of the elec- 
tric field using Eq. Comparing the injection rates 
given in Eqs. @ and ([3]) with the one-photon indirect 
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optical injection rates^ we find that these two formu- 
las differ only in the transition matrix elements given in 
Eq. (j3|). Therefore they show similar temperature depen- 
dence, which is mainly determined by the phonon num- 
ber, and similar contributions from each valence band, 
which is mainly determined by the joint density of states. 



To better understand these results, we first consider 
the properties of transitions around the band edges. 
Then we turn to the injection rates of carriers and spins, 
and the DSP under a~ light propagating along two dif- 
ferent directions. 



_ 0.4 
i 

> 
7 

to 

2s 0.2 
i — i 



i 1 1 i 1 1 i 1 1 i 1 y 

Cxxxx / 

SZ / 
txxyy / 

- ?Z / 

tzzxx / 
■ SZ / 


1 , 1 1 1 , 1 1 1 1 

c xyxy 

SZ 

txzxz 
- SZ 

. tT z 


/ y s 

/ 7 / 






,■// 





0.3 0.6 0.9 0.3 0.6 0.9 

2hw - E ig (eV) 

FIG. 1: (Color online) Spectra of £ a z bcd at 4 K (thick black 
curves) and 300 K (thin red curves). 
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FIG. 2: (Color online) Spectra of llm[C, f z abcd ] at 4 K (thick 
black curves) and 300 K (thin red curves). 

In a previous paper j- n we discussed in detail the photon 
energy, temperature, and phonon branch dependence of 
the total carrier injection coefficients f 111 , £ xx vv^ anc l 
{:xyxy _ Here the £2 in Fig. [T] show similar proper- 
ties: For excess photon energies 2hui — Ei g of interest, 
^zjixx g rs |. mcreases w ith increasing photon energy, and 
then slightly decreases; all the other components increase 
monotonically. In contrast, all the components of C,^ abcd , 
given in Fig. [2j show a complicated photon energy de- 
pendence. All injection rates at 300 K are larger than 
those at 4 K due to the larger phonon populations. 



1. Transitions near band edges 

One can try to simplify the description of the indirect 
two-photon injection around the band edges using the 
high symmetry at the band edge. We symmetrize the 
indirect two-photon injection rates as 
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Here A%J: kcVkvX — Hp cj •^-e(p c ,ik c )v{p e , I k v )x/-^' c > I 'i 
are the Af c ,i symmetry operations in d v that keep the 
I th valley unchanged, while P v are the Af v symmetry op- 
erations in Oh] and X)agt indicates summation over all 
modes in the r th branch. Around the band edge, it is a 
good approximation to take the mediated phonon energy 
Ml(k -kv)\ an( i the phonon number N( kc _ k ^\ to be con- 
stant and equal to their band edge values hfl k ox = hfl^ 
and N k o.\ — N® , respectively. Then the injection rates 
are approximately 
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Here the symmetrized expression Af?°j* vh T in Eq. (fTTj) 
avoids the ambiguity in calculating the band edge values 
of A a ^ d vk A , which is induced by the degeneracy of the 
heavy and light hole bands at the T points. This can be 
clearly shown by rewriting W^ kl , x = {ck c \W^ vkvX + 
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keeping the intermediate states appearing in Eq. (U) im- 
plicit. Then similar to the corresponding results for one- 
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photon indirect optical transition^ we have 

Jabcd Jabcd ^ \ Jabcd 

^/;cfcOHHfcgr — -^cfcOLHfcOr — 2 • A I;ck°v' k°T > 

i/=LH.HH 

which is unambiguous for any choice of heavy and light 
hole state at the valence band edge. We analyze the 
nonzero matrix elements of A^ k o vk o T using the symme- 
tries of the crystal. 

For a given symmetry operation, the transformation of 
W ab is determined by H^ p , v a , and v b ; a direct symmetry 
analysis for W^?o sfe o A is possible with the electron state 
|cfc") and the hole state \vk°). However, because of very 
weak spin orbit coupling in silicon, this process can be 
greatly simplified by dropping spin orbit coupling terms 
in H^, v a , and v b , and thus in vk y Without spin- 
orbit coupling, the valence states at T are chosen with 
the symmetry properties of {X = yz,y = zx,Z = xy}; 
the phonon states are chosen with the symmetry {x, y} 



for TA/TO branch, {z} for LA branch, and {x 2 - y 2 } 
for LO branch; without losing generality, the conduction 
band edge states are taken to lie in the Z valley, which 
has the symmetry of {z}. All matrix elements are listed 
in Table [T] In total there are fifteen nonzero quantities 
for the band edge values. From the table, selection rules 
depend strongly on phonon states. 

With spin orbit coupling, the valence bands are split 
into HH (||,±|)), LH (|§,±±)) and SO (|±,±±)) bands, 
and the conduction bands are two-fold spin degener- 
ate bands \z f) and \z I). The indirect optical ma- 
trix elements in these states can be easily obtained by 
linear combination of the terms in Table U and the 
band edge transition probabilities can be identified by 
Af£ T = Af k f vk o x , with k° c being the band edge wave 
vector in the Z valley. Similar to the corresponding term 
for one-photon absorption,— A'^'. c c f, T has the following 
properties: i) ^-cvt are the same for v = HH, LH, and 
cjo- iil V ?f abcd — n nnH rf abcd — ?f abcd 

Bu > "J 2^v t>Z;cvT — U ana SZ;cHHr - Sz ;c LHt- 
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TABLE I: Band edge value of W^o vk o x - Here M l] is a matrix with matrix elements [M lJ ]ki — (l — Sij)(Sik8ji+5aSjk) + 8ij5ik5u. 
There are totally fifteen parameters {Wi, i = 1, • • • , 10} for TA, LA, LO phonon branches and {IT/, i = 1, • • • , 5} for the TO 
phonon branch. 



We list -4^cHHr m Table |TT] for carrier injection and 
Table IIIII for spin injection. Generally, these nonzero 
transition probabilities can be used in Eq. (jlip to ap- 
proximate the Aj b ^ k vk T around the band edge values, 
which results in a simple formula 

2"7T 

A ab lt± « -j-J cv (hu)N? ± A a j^ T , (14) 

the analog of which is widely used in the qualitative anal- 
ysis of one-photon direct and indirect injection even for 
injection away from the band edge. Here J cv (huj) is the 
joint density of states, J cv {hw) = J2k c &i-k v ^( £ ck c -£vk^± 
ftQ,® — fvjj). In Fig. EJa), we give the local properties of 
^"c"(HH))t„T around band edges (fc°, fe°); its rapid varia- 
tion away from the band edge shows that the simple for- 
mula (fT4]) may fail. Garcia and Kalyanaramara2£ found 
that the corresponding formula for two photon absorp- 
tion should be replaced by 
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TABLE II: Independent nonzero components of £z?ci 



Here the two-photon absorption coefficient (3 is re- 
lated to our calculated quantity £ xxxx by f3 = 
2hui^ xxxx I (2nRceo) 2 , jib. is the refractive index, c is the 
speed of light, eo is the vacuum permittivity, F n (x,y) — 
(x — y — l) 2+n /x 5 , and = fl k a\ is the frequency of 
phonons mediated in the band edge transitions. Ac- 
cording to the parity difference between the band edge 
hole and electron states, the transitions are divided into 
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FIG. 3: (color online), (a) Values of £z Z cfc c (HH)fc„(TA) a ^ 
2&j = Eig along different directions: (Dashed black curve) 
(k° + kz,k°), (Solid red curve) (fe° + ky,k°), and (Dotted 
blue curve) (fe°,fc° + kx). (b) photon energy dependence 
of (2fc) 6 Cz Z cfe°(HH)fe°T f° r different phonon branches. Here 
a — 5.43lA is the lattice constant. 



allowed-allowed (a — a), allowed- forbidden (a — f), and 
forbidden-forbidden (/ — /) processes, which correspond 
to the n = 0, 1, and 2 terms in Eq. (flSl . respectively; such 
a classification is based on whether the band edge values 
of the matrix elements of v a and v b in Eq. (j4} are zero 
(forbidden) or nonzero (allowed) for different parity of 
the intermediate states. In deriving Eq. (|15p . the depen- 
dence £ckcvkv\ (fi^) 6 must be used. Dinu-5, argued 
instead that ^k^vk v x ^ (^-")~ 5 for some processes. Here 
we can numerically study this dependence, and the result 
is plotted in Fig.[3][b); it shows that w -5 and u~ 6 depen- 



dences are both important, at least for the TA phonon 
branch. 

The simple formula (|14l) only corresponds to the a — a 
process. In obtaining results for the other two processes, 
the first and second derivatives with respect to k c and 
k v of Wf^ ck vh x are necessary. This results in a more 
complicated symmetry analysis that we do not consider 
here. 

2. Injection for a ~ light propagating along (001) and (111) 

For <j~ light propagating along the directions (001) and 

(111), the electric field E^ can be written respectively 
as 



p(ooi) _ ^ *y p 

E " ~7T Eo ' 



E (in) = + (yg - i)y - (V3 + i)z 



2V3 



E , (16) 



where (k) denotes (001) or (111). The injection rates of 
carriers and spins then are 



n I ., cv x± = Ci k Lx±\ E o\ 4 » 

of _ /-f\(k) I 77. |4 

b I;cvX± - l =I;cvX±\ £j 0\ ■ 



(17) 



Here and are the injection coefficients of car- 

riers and spins, respectively, in the I th valley. They can 
be expressed by the nonzero components of £,z h -cvX± anc ^ 
Cz a cvX± defined in Eqs. (HJ and ((6|). The resulting ex- 
pressions are listed in Table IIV1 The DSP is defined as 
DSP / = S f /(hh/2). For (001) light, the injected spin 
in each valley and the total spin are all parallel to the 
light propagation direction, i.e., the (001) direction. The 
carrier and spin injection rates show valley anisotropy be- 
tween the Z valley and the X/Y valleys. For (111) light, 
the injected carriers are the same for every valley, and the 
total spin polarization is still along the direction of the 
electric field, i.e., the (111) direction. But the injected 
spins in each valley have different spin polarization: The 
two transverse directions in each valley have the same 
injection rates, which are different from the longitudinal 
direction of the valley. 



3. Carrier injection under o light propagating along (001) 
and (111) 

We plot photon energy dependence of the total carrier 
injection coefficients for (001) and (111) light at 4 K and 
300 K in Fig. [4] (a). As we found earlier in a prelimi- 
nary study^ the injection coefficients increase rapidly 
with increasing temperature. The injection for (111) 



light is larger than that for (001) light, demonstrating the 
anisotropy of the injection on light propagating direction. 
In agreement with Hutchings and Wherrett's notation^! 
for direct gap two-photon injection, this anisotropy can 
be characterized by two parameters, the anisotropy a and 



7 



(fe) 


/ 


t(fc> 


t x;{k) 
Si 


si 


SI 


<ooi) 


X 


1 / szzzz 1 1 t-zzxx . 
4"ltz + ?Z ) — 2SZ + S"Z 








MCz M + Cz ] 


Y 


e(00l) 

SX 








*z;(001) 


z 


1 fxrai 1 rxxyy . c xyxy 

-?z ~ -?z T SZ 








2Im[C^ !,M ] 


Total 


1 £xxxx 1 £xxyy _■_ £xyxy 
2"S ~~ 2"S "T S 


o 







(111) 


X 


d in> 




^;(lll> 


^;{lll) 


Y 


sz 


sz 


SZ 


sz 


Z 


1 -^Z — <^z — ^z 
9 +4 ( € »W +2f g aaI *) 


2 t r ^xyzyy . nA^^z^i/ i >xzzyzi 

57f MCz + 2 Cz + Cz ] 


r x;(lll) 

sz 


4 t \ j-zxvxx i /-Z1VZXZ1 

s^sMCz + Cz" ] 


Total 


1 ^aessz _ ^xxyy j_ ^xyxy^ 






^;{111) 



TABLE IV: The carrier indirect two-photon injection coefficients £j and the spin indirect two-photon injection coefficients 
(-a,(k) ^ e jth va ]\ e y f or G - light propagating along directions (k) — (001) and (111). 
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FIG. 4: (Color online) (a) Spectra of total carrier injection 
rates £ at 4 K and 300 K for <r~ light propagating along 
the (001) and (111) directions, (b) Anisotropy a and linear- 
circular dichroism 8 at 4 K (thick black curves) and 300 K 
(thin red curves). 
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FIG. 5: (Color online) Spectra of carrier injection rates ^j 001 ^ 
in the I = Z and X valley at 4 (thick black curves) and 300 K 
(thin red curves). 



the linear-circular dichroism S, which are given as 

£xxxx <2£:xyxy pxxyy 



Cxxxx 

2£xyxy _|_ frxxyy 



2£ a 



(18) 



In the isotropic limit, a — and S — ^ xx vv / ^ xxxx ^ 
We plot cr and 5 in Fig. @] (b). Note that the anisotropy 
shows a much stronger temperature dependence than the 
linear-circular dichroism. In contrast to a in direct gap 
two-photon injection, which clearly shows the onset of the 
transition from the spin split-off band to the conduction 
band by the presence of a cusp, it is hard to identity 
the contribution from spin split-off band in indirect gap 
injection. This is because the energy dependence at the 
onset of indirect absorption is oc (hw — Ei g ) 2 , given by 
the a — a process, instead of oc (hu — Eg) 1 / 2 for direct 
absorption. 

Now we turn to the carrier injection into each val- 
ley. For (111) light, all valleys are equivalent, and the 
injection coefficient in each valley is 1/6 of the total. 
There is no valley anisotropy in this case. For (001) 



light, the valleys can be divided into two sets: {Z, Z} 
and {X, X, Y, Y}, and the injection is the same for all 
valleys within each set. We plot the spectra of injec- 
tion rates £j 001 ^ in the I = Z and X valleys at 4 K and 
300 K in Fig. [5J The spectrum in each valley has a shape 
similar to the total, and the injection in the Z valley is 
larger than that in the X valley. The valley anisotropy 
arises because of the anisotropic effective electron mass 
in the conduction bands, which leads to different matrix 
elements appearing in (j4]) for the different Cartesian com- 
ponents of velocity. For the Z valley, the effective mass 
along the z direction is heavier than that along the x/y 
directions^ which results in a smaller z component of 
the interband velocity matrix elements.— From Table ITVl 
we find that the z-components of the electron and hole 
velocity only appears in the injection rates in the X/Y 
valleys, and results in their smaller injection rates. 

Figure [5] gives the phonon-resolved spectra in the X 
valley for (111) light. Similar to our previous results^ 
we find here that the LA phonon-assisted process gives 
the smallest contribution, while the TA and TO phonon- 
assisted processes dominate: At low temperature, the 
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FIG. 6: (Color online) Phonon resolved spectra of carrier in- 
jection rates f at 4 K (thick black curves) and 300 K (thin 
red curves). 



TA phonon-assisted process dominates at low photon en- 
ergy, and the TO phonon-assisted process dominates at 
high photon energy; with temperature increasing, the TA 
phonon-assisted process becomes more and more impor- 
tant due to the small TA phonon energy, and dominates 
for photon energy less than Ei g at 300 K. The phonon- 
resolved injection rates in each valley for (001) light show 
similar behavior. 



4- Spin injection under a light propagating along (001) 
and (111) 

In Fig.[7]we show the spectra of the spin injection rates 
£^ and the DSP-* for o~ light propagating along (001) 
and (111) directions at 4 K and 300 K. The total spin 
polarizations are all parallel to light propagation direc- 
tion. When photon energy is higher than the injection 
edge, which is E ig + Ml^ A at 4 K or E ig — hQ^ Q at 300 K, 
the spin injection rates first increase with photon energy 
from zero to maximum values, then decrease, and then 
change direction at high photon energies. This is different 
than the behavior of indirect one-photon spin injection, 
in which the spin injection rates always increase with 
photon energy. The difference can be attributed to the 
complicated transition amplitude W in Eq. The fine 
structure of the injection rates around the band edge are 
clearly shown by the DSP spectra in Fig.Efb). The DSP 
depends strongly on the laser propagation direction and 
the temperature. For (001) light, the maximum DSP can 
reach about 20% at 4 K and -10% at 300 K; for (111) 
light, the maximum DSP is only 6% at 4 K but 20% at 
300 K. Around the injection edges, the DSP show more 
detailed structures at 300 K than 4 K. In Fig. EJb), we 
label the injection edge for phonon branches by dotted 
vertical lines: hui, hu>2, and hu>3 for the TO, LO, and TA 
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FIG. 7: (Color online) Spectra of (a) spin injection rates Q k ' 
and (b) DSP (fe> at 4 K (black curves) and 300 K (red curves). 
(Solid curves) E || (001), (dashed curves) E || (111). The spin 
polarization direction is parallel to E. The labeled energies 



are 2hw\-z 



E ig 
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phonon absorption process, respectively; Hu}' 3 identifies 
the TA phonon emission process. As in the correspond- 
ing results for indirect one-photon spin injection, the fine 
structures arising here come from the contributions of 
different phonon branches. 

To better understand these fine structures, we plot the 
spin injection in each valley and the contribution from 
each phonon branches for (001) light in Fig. [5J Fig- 
ure [5Ja) gives the spin injection rates in the X and Z 
valleys at 4 K and 300 K, in which the valley anisotropy 
is prominent. Again, due to the anisotropy in electron 
velocity, the injection rate in the Z valley is much larger 
than that in the X valley. Figure H](d) gives the detailed 
structure of DSP around the band edge. The maximum 
DSP is about 36% at 4 K and -20% at 300 K in the X 
valley, and about 10% in the Z valley for both tempera- 
tures. However, the spin injection rates are very close in 
these two valleys near the injection edge in Fig. [T^a), so 
the difference between these maximum values can only 
come from the difference of the carrier injection rates, 
which are much smaller in the X valley than that in the 
Z valley (see Fig. [5J. 

The phonon-resolved spin injection rates in the X and 
Z valleys are plotted in Figs. [8] (b-c) at 4 K . In the X 
valley, the TA phonon branch dominates at low photon 
energy and the TO phonon branch dominates at high 
photon energy. In the Z valley, the TA and LO phonon 
branches have similar contributions and dominate at low 
photon energy. Near the band edge, the spins injected 
from the TO and TA phonon-assisted processes have op- 
posite spin polarization direction in the X valley, but 
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FIG. 8: (Color online) Spectra of (a) spin injection rates ^j 001 ^ 
in the I — X valley (solid curves) and the I = Z valley 
(dashed curve) at 4 K (thick black curves) and 300 K (thin red 
curves), (b)-(c) phonon branch resolved spin injection rates in 
the X and Z valleys at 4 K, and (d-f) The corresponding DSP 
of (a-c).The labeled energies are 2fiwi_ 2 = and n£ A , 

respectively. 



same in the Z valley. Figures [8] (e-f) give the corre- 
sponding DSP. Almost all processes contribute nonzero 
DSP. In the Z valley, the spin injection rates are given 
by Im[£jj^ a:a! ] as shown in Table HVl At the band edges, 
the carrier and spin injection amplitude of the LA and 
LO phonon-assisted processes are all zero, which means 
these a — a processes inject no carriers. As one moves 
away from the band edges, carriers and spins can be in- 
jected by a — / and / — / processes, which results in a 
nonzero DSP. In the X valley, the spin injection rates 
are given by Im[(^" !,z + C^.^. vv \. From the results in Ta- 
ble [Ell we see that the LA phonon-assisted process gives 
zero spin injection amplitude at band edge, but its DSP 
is not zero because the a — f and f — f processes domi- 
nate over the a — a process. Similar results also exist in 
two-photon direct injection. 50 

In Figs. [8] (e-f), we plot only the DSP for the phonon 
emission processes, the corresponding injection edges are 
given by Huj^. From the calculation of the indirect one- 
photon injection, we know that the DSP induced by the 
phonon absorption process and the phonon emission pro- 
cess have a similar shape, but the injection edge shifts 
from Su^ to hwi. At 4 K, the injection edge is domi- 
nated by the TA phonon emission process (which begins 
at Huj' 3 ), and it is dominated by the TO/LO phonon ab- 
sorption process (begins at hu>i) at 300 K, then follows 



by the LA phonon absorption process at ftioi and the TA 
phonon at Hui^. Therefore, the co-action of the TO and 
LO phonon absorption processes gives the negative DSP, 
and results in the sharp increase between the photon en- 
ergies huii and Hui2 in Fig.[8](c), then the LA/TA phonon 
absorption processes give positive DSP, so the total DSP 
decreases sharply after hio^. 

Figure [S] gives the details of the spin injection for the 
o~ light propagating along (111) direction. The analysis 
is similar to the (001) case. 
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FIG. 9: (Color online) Spectra of (a) spin injection rates 

C# and (d) DSP^ 111 ^ for the a — z (solid curves) and 
x, y (the dashed curves) spin components at 4 K (thick black 
curves) and 300 K (thin red curves). Phonon resolved (b- 

c) spin injection rates f^'.^ 111 ^ and (e-f) DSP^ 111 ^ f° r the z 
(black thick curves) and x, y (blue thin curves) spin compo- 
nents at 4 K. 



III. TWO-COLOR CHARGE AND SPIN 
CURRENT INJECTION 



Now we study the motions of optically injected carriers 
and spins. While a single color light source can inject 
net current into a particular valley ) 51 ' 52 due to the Oh 
symmetry there is no net charge or spin current injection 
from either either one-photon or two-photon absorption 
of a single color light source. We calculate 1 + 2 injection 
effects here, and only consider the total charge and spin 
current induced. For a two-color optical field E(t) — 
-EL,e~ iwt + i?2ue~ j2ajt + c.c, the carrier density injection 
rate is 



Z™ezk(KKY , (is) 
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where £f & are one-photon indirect injection coefficients 39 
and £ abcd are the two-photon injection coefficients studied 
in the previous sections. The interference between u> and 
2lo beams injects charge and spin currents with injection 
rates 



j d 



k fd 



(h) 



vit^EUElEZy +c.c 



dfabc EUEtEZT +c.c. 



(20) 



The injection coefficients r] d ^ and V d ^ c are written in 



the form B 



abc 



Ea±«a± ^th 



be 



na, 
D A± 



vjabc 

^ck c vk v X 



ck c ,vk v 



= » 2J (c / fec|(«'fet;|'B|wfe t ,)|cfe c 

x T,f fecCfeuA (2 W )[M/|, c fccC ,^ A ( W )]* . 



(21) 
(22) 



Here T° fc - k x is the one-photon indirect optical transi- 
tion amplitude^ 9 - By taking B as J d = —ev d , Jfi — ev d , 
K[ d = ~i(v e d Si + Sfv d e ), and k{ d = ^(y h d s{ + S{v d ), 
we obtain the injection rates for electron and hole charge 
and spin currents, with rj dabc = rj dacb and ^ ahc = ^ rf / Qcb . 
For diamond structure crystals, the nonzero components 
are 
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(23) 



(24) 



resolved tensors ^(/if-Ai and M e 
symmetry properties as the total injection tensor r\ 



All other nonzero components can be obtained by cyclic 
permutations of the Cartesian indices. The phonon- 

"t{h);\± and M*$Lfc Share the Same 

dabc 
e(h) 

and Mg|^j fcc , respectively. Using time-reversal symmetry, 
in the independent particle approximation we adopt here 
all r] dabc are pure imaginary numbers, and all ^, d f abc are 
real numbers. We show the calculated spectra of each 
component of the charge current in Fig. 1101 and of the 
spin current in Fig. 111! The current injection coefficients 
^e(h) and ^e(h) C nave the same symmetry properties as 
that of the two-color direct current injection across the 
direct gap of germanium^ 
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FIG. 10: (Color online) Spectra of Im[ri fal " : ] for electron (a) 
and hole (b) at 4 K (thick black curves) and 300 K (thin red 
curves) . 
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FIG. 11: (Color online) Spectra of /if /a6c for electron (a) and 
hole (b) at 4 K (thick black curves) and 300 K (thin red 
curves) 



From the calculation, both the charge and the spin 
currents for injected electrons are larger than that for in- 
jected holes. One contributing factor is that the electron 
moves faster than the hole due to the smaller effective 
mass. But for the spin current, another factor is that the 
average spin expectation value over the HH and LH band 
is smaller than that in the conduction bands. 

We consider the indirect current injection coefficients 
under the two-color laser beams propagating along the 
z direction with the electric field components taken as 
E u = Eue^-eu and E 2uJ = E 2ui e i ^e 2ui . Here E u and 
E 2 uj are real and positive field amplitudes, and e 2uJ 
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are polarization vectors, and 4>u and 02 W are their phases, 
and A<p = 2<j) u — (p 2u > is the relative-phase parameter that 
is used to control the current. In the following, we give 
the current injection for different configurations of the 
laser beams. 



A. Co-circularly polarized beams 



For two circularly polarized beams propagating along 



the z direction, the electric fields are e 2u 



r^2 



and 



e u = a si with Si — ± identifying for the handedness and 
<t s = (x + isy) I y/2. The indirect gap current injection 
coefficients are 
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iui z a m b (25) 



with m Sl = s\x sin Acj> + y cos Acp. Both the direction 
of the charge and spin currents and the polarization of 
the spin current can be controlled by the relative-phase 
parameter A<p and the light polarization s$. The charge 
current flows only in the x — y plane, and the calculated 
^e(h) IS negligible. For the opposite circularly-polarized 
beams, the x component of the charge current is kept 
unchanged, but the y component reverses. The spin cur- 
rent flows in the x — y plane with spin polarization along 
x-axis, or flows along z-direction with spin polarization 
along x/y direction. 



B. Cross-linearly polarized beams 

For two z propagating cross-linearly polarized beams, 
along the x direction and E 2u j along the y direction, 
the injection current rates are given as 



jy-ab _ 
K e(h) ~ 



2lm[ n ™™]ElE iu ysmA<t> , 

2(m; 



(26) 



zxyxx ~a 
e(h) 



b zxyzz * a 
x -M e(h ) 



.1 *z b )ElE 2u cosA</> 



In this scenario, the charge current and the spin current 
are injected with it/2 phase difference. Therefore, by tun- 
ing the relative-phase parameter A<f>, a pure charge cur- 
rent or pure spin current can be injected. The charge cur- 
rent flows along the second harmonic polarization axis, 
and its amplitude is determined by rj xxyy , which is zero 
under the parabolic band approximation. In our full 
band structure calculation, it is nonzero due to the band 
warping, but very small compared to other tensor compo- 
nents. The spin current has two components, one involv- 
ing flow along the x direction with the z spin polarization, 
and the other involving flow along the z direction with 
the x spin polarization. 



C. Co-linearly polarized beams 

For two z propagating beams, both polarized along the 
x direction, the injection current rates are given as 



Je(h) 

Tsab 
K e(h) 



2lm[r, x e fff]ElE 2uXS mA<j> , 
2^\y a z h - z a y b )ElE 2u cos A0 . (27) 



This scenario also gives the phase difference between the 
charge current and the spin current as 7r/2, so as for cross- 
linearly polarized beams pure charge current injection or 
the pure spin current injection can also be realized by 
choosing a suitable relative-phase parameter A<j>. Our 
results give the relative-phase parameter dependence of 
the injected current as sin A(f>, which is in good agreement 
with the experimental results^ around zero probe delay. 
To understand the indirect current injection better, we 
compare the indirect current injection with the direct 
one. Because of the lack of the direct gap injection in 
silicon in the literature, our results are compared with 
the direct current injection in bulk germanium.^ For the 
charge currents injected across the indirect gap in silicon, 
the electron and hole currents have opposite directions at 
high photon energies, but they can be the same at low 
energies; for charge currents injected across the direct 
gap in germanium, they always have the same directions. 
For the spin current, the injected spin current is not so 
small compared to other components, especially at 300 K, 
while they are ignorable small in the direct gap current 
injection in germanium because of the complete lack of 
the helicity of the incident light. 

In this configuration, a good characterization of the 
charge current is the swarm velocity, which is defined 
as the average velocity per injected carriers forming this 
current, v x = J x /e s n, with h taken from Eq. (fT9|) . Here 
e s = —e is used for electrons and e s = e for holes. 
When A<j) is a multiple of 717 '2 and the indirect one- 
photon charge injection rate equals the indirect two- 
photon charge injection rate, the maximum swarm ve- 
locity is 



lm[r] xxxx (oj)] 
e s ^/^ xx (2uj^ xxxx (u J ) 



(28) 



We show the maximum swarm velocity in Fig. [12] for 
the injected electrons and holes at 4 K and 300 K. The 
behavior of the swarm velocity can be divided into two 
regions: (i) for photon energies in 2Huj — Ei g ^ 0.25 eV, 
the maximum swarm velocities are along the x direction 
for injected electrons and holes, and both magnitudes in- 
crease with increasing photon energy. Compared to the 
maximum swarm velocities in bulk germanium, which is 
in the order of 10 3 km/s, the velocity here is about one 
order of magnitude smaller due to the larger conduction 
band effective mass in silicon, (ii) for photon energies in 
2hw — Ei g < 0.2 eV, the swarm velocities show fine struc- 
tures. In particular, all currents experience directional 
changes except the electron swarm velocity at 4 K. Anal- 
ogous to indirect one- and two-photon charge and spin 
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FIG. 12: (Color online) Maximal swarm velocity max for 
the injected electrons (solid curves) and holes (dashed curves) 
at 4 K (thick black curves) and 300 K (thin red curves). 
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FIG. 13: (Color online) Phonon-resolved Maximal swarm ve- 
locity for the injected (a) electrons and (b) holes at 4 K (thick 
black curves) and 300 K (thin red curves). 



injection, these fine structures are induced by the dif- 
ferent phonon branches, which are clearly shown in the 
phonon-resolved maximal swarm velocity in Fig. 1131 At 
300 K, the injection edge is given by the TO/LO phonon 
absorption processes, both of which give negative veloci- 
ties for injected electrons and positive velocities for holes 
around the injection edge. At 4 K, the injection edge is 
given by the TA phonon emission process, which direc- 
tion is opposite to the band edge current at 300 K. There- 
fore, the sign change of the injected current is induced 
by the contributions from different phonon branches. At 
high photon energies, the injected velocities are almost 
independent of the temperature. This is similar to the 



temperature dependence of the DSP of the one-photon 
indirect injection^: The only temperature dependence 
in the injection rates lies in the phonon number, which is 
the same for the denominator and numerator in Eq. (|28[) 
for a given phonon branch. For high photon energies, an 
average phonon number can be used as a good approx- 
imation, and the swarm velocities, given by the ratio in 
Eq. (|2"5)l . are almost temperature independent. 



IV. CONCLUSION 

In conclusion, we have performed a full band structure 
calculation of two-photon indirect carrier and spin injec- 
tion, and two-color indirect current injection, in bulk sil- 
icon. We presented the spectral dependence for all com- 
ponents of the response tensors at 4 K and 300 K, with 
which the injection under any laser beams can be ex- 
tracted. All injection rates increase with increasing tem- 
perature due to strong electron phonon interaction at 
high temperature. We discussed in detail the injection 
under different polarized light beams. 

For two-photon indirect optical carrier and spin injec- 
tions, we considered the injection under a~ light propa- 
gating along (001) and (111) directions. For (001) light, 
the injection rates in the X and Y valleys are the same, 
but different from that in the Z valley; for the (111) light, 
the injections into all valleys are equivalent. For car- 
rier injection, injections for these two light propagation 
axes differ slightly. The calculated injection anisotropy 
and the linear-circular dichroism characterize the non- 
parabolic band effect in the full band structure calcu- 
lation. For the (001) light, the injection in the Z val- 
ley are much larger than that in the X/Y valleys, and 
give the valley anisotropy, which is induced by the veloc- 
ity anisotropy in the conduction band. At 4 K, the TA 
phonon-assisted process dominates at low photon ener- 
gies, and the TO phonon-assisted process dominates at 
high photon energies. At 300 K, the TA phonon-assisted 
process dominates for all photon energies. 

For two-photon indirect gap spin injection, the total 
injected spins orient parallel to the light propagation di- 
rection for the two directions considered. The spin in- 
jection rates increase from the injection edge to a maxi- 
mum value with the photon energy increasing, and then 
decrease. The DSP strongly depends on the tempera- 
ture around the injection edge. For the (001) light, the 
injected spins in each valley are still along the z direc- 
tion, but the spin injection rates in the X and Z valleys 
are different. The maximum DSP of total spins is 20% 
at 4 K and -10% at 300 K; the DSP can reach about 
40% at 4 K and -20% at 300 K in the X valley, and 
both become 10% in the Z valleys. For (111) light, the 
spins in each valley orient to a direction different from the 
light propagating direction. In the Z valley, the x and 
y (transverse) components have the same injection rates, 
which are different from the z (longitudinal) component. 
The maximum DSPs of the total spins are 6% at 4 K and 



13 



20% at 300 K, while the one of the z component spin in 
the Z valley is about 5% at 4 K and 20% at 300 K, and 
becomes 25% and —8% for the x or y components. All 
these features are induced by the interplay of different 
phonon branch-assisted processes. 

For light propagating along the (001) direction, the 
injected carriers or spins break the symmetry between 
the X and Z valleys. Such a valley anisotropy of injected 
carriers could be probed experimentally, for example, in 
a pump-probe scenario where the probe beam propagates 
cither parallel or perpendicular to the pump beam. 

For the coherent control, we calculated two-color in- 
direct charge and spin current injection under three dif- 
ferent polarization configuration of the two-color beams 
propagating along the z direction. For the co-circularly 
polarized beams, the direction of the injected charge cur- 
rent is in the x — y plane; the spin current flows in the 
x — y plane with a z oriented spin polarization, or flows 
along the z direction with the spin orientation in the 
x — y plane; the current direction or the spin polariza- 
tion in the x — y plane can be controlled by a relative- 
phase parameter. For the co-linearly polarized beams 
and the cross-linearly polarized beams, the directions of 
the charge current, the spin current, and the spin polar- 
ization are orthogonal to each other. In these two cases, 
a pure spin current or a pure charge current can be ob- 
tained by choosing a suitable relative-phase parameter. 
We calculated the maximum swarm velocity for charge 
current as a function of photon energy, and found that 
the maximum swarm velocities undergo a sign change 
near the band edge, which is induced by the contribu- 
tions from different phonon branches. 
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Appendix A: Dependence of the injection rates on 
the light propagating direction 

For a circularly polarized laser pulse propagating along 
direction fi\, the electric field can be expressed as 



El E . „ . 

E u = — (n 2 + tsn 3 ) 
v2 
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here s = ±1 identifies the helicity. In the Z valley the 
carrier injection rates can be written as 
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with fi\ = z x fi2 and 
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gi (e, (j>) = (shr 6»sin40)n 2 + [4 + sin 6>(cos^ 

. ' (A " 5) 

For an arbitary propagating direction (9, </>), the direction 
of the spin polarization is not always along hi. In this 
case, we define DSP as the ratio of the magnitude of the 
spin injection rate and the carrier injection rate P z = 
\S z \/(n z h/2). We plot in Fig.[TJ]the (0, ^-dependence 
of the P z at the edge of the each phonon emission process 
in the Z valley, which shows strong anisotropy of the 
light propagating direction. The maximum P z can reach 
45% at (0, <t>) w (7r/2,0.27r) for the TA phonon emission 
process, 20% at rs (7r/2,0.l7r) for LA phonon, 45% at 
« vr/4 for LO phonon, and 13% at w 0.2tt for TO 
phonon. 

The total carrier injection rates are 



EU 



4 /-xxxx 



l 



■ sin 2 0(sin 2 <f> cos 2 <j> sin 2 + cos 2 0) 



(A6) 



with S and a defined in Eqs. (|18l) . The total spin injection 
rates are 

S = £#sj - 2lm[C xyxx ]h 1 

+ jsm9g 2 (9,cf>)Im[( zxyxx - 2( xxzx y}} , (A7) 



with 



92(0,0) = sin 0(1 + 7 cos 2 - sin 2 cos 40)n.i 

— (sin 2 0sin40)n2 

- cos 0(3 - 7 cos 2 + sin 2 cos 4^)n 3 .(A8) 

n 3 — nixri2 , p rom above expressions, the two-photon carrier and spin 
injection rates show strong anisotropy for the light prop- 
(A2) agating direction. For the spin injection, the direction 
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6 X 7T 

FIG. 14: (Color online) Light propagating dependence of the 
P z at the edge of the each phonon emission process in the Z 
valley, (a) TA phonon, (b) LA phonon, (c) LO phonon, (d) 
TO phonon. 



of the injected spin polarization usually differs from the 
light propagating direction, but it reverses as the light 
helicity changes. 

For the linear polarized laser pulse, we also found that 
the total carrier injection rates strongly depend on the 
polarization direction. For the electric field 

= E (n 2 sin a + n 3 cos a) , (A9) 

with a for the polarization direction, the total carrier 
injection rates are given by 

n = £ 4 £ OT {1 - [1 - f(9, & 0)]a} , (A10) 

with 

/(#, 0, /3) — (cos <f> sin a — cos 6 sin <f> cos a) 4 
+ (sin <j) sin a + cos 6 cos cf> cos a) 4 
+ sin 4 9 cos 4 a. (All) 
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